Effect of heating and cooling of photogenerated electron-hole plasma in optically 

pumped graphene on population inversion. 
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We study the characteristics of photogenerated electron-hole plasma in optically pumped graphene 
layers at elevated (room) temperatures when the interband and intraband processes of emission and 
absorption of optical phonons play a crucial role. The electron-hole plasma heating and cooling 
as well as the effect of nonequilibrium optical phonons are taken into account. The dependences 
of the quasi-Fermi energy and effective temperature of optically pumped graphene layers on the 
intensity of pumping radiation are calculated. The variation of the frequency dependences dynamic 
conductivity with increasing pumping intensity as well as the conditions when this conductivity 
becomes negative in a certain range of frequencies are considered. The effects under consideration 
can markedly influence the achievement of the negative dynamic conductivity in optically pumped 
graphene layers associated with the population inversion and, hence, the realization graphene-based 
terahertz and infrared lasers operating at room temperatures. 
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I. INTRODUCTION 



The eapless energy spectrum of graphene layers 
(GLs) [i| provides an opportunity to create devices uti- 
lizing the interband transitions in a wide spectral range 
of radiation. As shown (see also Refs. 0-111), under 
optical or electrical pumping, the interband population 
inversion can be realized in GLs. Due to fairly high quan- 
tum efficiency of the interband transitions in GL 0| 
and particularly in multiple-GL structures [1, Q , the in- 
terband contribution to the real part of the dynamic con- 
ductivity can surpass that associated with the intraband 
(Drude) absorption [21] , so that the GL net dynamic con- 
ductivity is negative in a certain range of signal frequen- 
cies. This opens up the prospects of achieving of las- 
ing of the terahertz (THz) and infrared (IR) radiation in 
different structures on the base of single- and multiple- 
GLs [Tol - [l^ . In optically pumped GLs, a marked por- 
tion of the absorbed optical energy can go to the elec- 
tron and hole heating. As demonstrated (sl. [l3l - [T5j . the 
electron and hole heating hinders the realization of pop- 
ulation inversion and negativity of dynamic conductivity. 
At elevated (room) temperatures, in contrast to the low 
temperature situation jl3| , the interband transitions as- 
sociated with optical phonons can be the dominant mech- 
anism of the recombination [l6j . The energy relaxation 
of electrons and holes in GLs can also be mainly due to 
interaction with optical phonons ^17l |. In such a case, 
the optical photons emitted by the photogenerated elec- 
trons and holes can accumulate in GL, so that the optical 
phonon system is heated as well. This can lead to the op- 
tical photons reabsorption accompanied with the intra- 
band electron and hole transitions. The deviation of the 



optical phonon system from equilibrium can also affect 
the properties of GL, as it occurs in carbon nanotubes 
(see for instance Refs. [l8l - l2^ ). in particular, affect the 
interband generation-recombination processes. Increase 
of the number of nonequilibrium phonons result in the 
increase of stimulated optical phonon emission and may 
lead to the generation of coherent phonons in graphene 
like it was discussed for a different systems [25-27]. Thus, 
the accounting for the optical phonon heating appears to 
be indispensable. However, it is important to stress that 
at certain values of the pumping radiation frequencies, 
cooling of the electron-hole plasma below the tempera- 
ture of the lattice is possible. 

In this paper, we study the effect of heating and cool- 
ing of the electron-hole plasma in optically pumped GLs 
at elevated temperatures, when the interaction with opti- 
cal phonons is the main mechanism of the recombination 
and energy relaxation, on its characteristics, particularly 
on its dynamic conductivity and in the THz and IR fre- 
quency ranges. The deviation of the optical phonon sys- 
tem from equilibrium (optical phonon heating) is taken 
into account The obtained characteristics are crucial for 
the realization of THz and IR lase rs Wi as well as differ- 
ent THz plasma wave devices [28l . [29l| on GL structures 
operating at room temperatures and can be used for the 
device optimization. 



II. MODEL AND THE PERTINENT 
EQUATIONS 

We consider the interband phogeneration of electrons 
and holes by optical radiation with the photon energy 



2 



Hft (optical pumping), where h is the reduced Planck 
constant. After the photogeneration of an electron and 
a hole with the energy exceeding the optical phonon 
energyfiwo ~ 200 meV, their behavior can correspond 
to the following tracks: 

(a) If the characteristic time, Tcc, of inter-carrier scat- 
tering (electron-electron, electron- hole, and hole-hole) is 
much larger than the characteristic time, tq of sponta- 
neous emission of optical phonon just photogenerated 
carriers can manage to emit cascades of optical phonons. 
The number of optical phonons K in such a cascade is 
determined by the ration n/ujQ. Due to the symmetry of 
the energy spectrum, the photogenerated electrons and 
holes emit equal number (K) of phonons. This cascade 
is followed by the electron and hole " fermisation" and re- 
combination i2] , so that the electron and hole distribution 
function becomes equal to / = {exp[(£ — ei?)/T] + 1}^^. 
Here e — vwP is the energy of electrons and holes with 
the momentum p, vw = 10® cm/s is the characteristic ve- 
locity of the GL energy spectrum, e^? is the quasi- Fermi 
energy, and T is the effective temperature (in the en- 
ergy units). The quasi Fermi energy of the electron- hole 
plasma in the situation in question ep and its effective 
temperature T are determined by the rates of photogen- 
eration ,and recombination, as well as the rates intraband 
and interband energy relaxation. These parameters are 
determined by the energy eg which is supplied to the 
electron-hole plasma in each act of the photogeneration. 
In the case in question, Sq = MIq = Ml — 2KhujQ. Gener- 
ally, sp and T are not equal to their equilibrium values, 
i.e., Ef 7^ and T ^ Tq; 

(b) If the inter-carrier scattering time is the shortest one 
(tcc ^ Tq), just photogenerated electrons and holes are 
immediately " fermilised" . In this case, all the energy of 
photogenerated carriers goes to the electron-hole plasma, 
so that So =^ hi}; 

(c) At relatively low energies of photons {hil < 2fia;o), the 
situation is similar to the previous case with Eq = hfl, but 
for rather arbitrary relationships between the character- 
istic times of inter-carrier and optical phonon scattering. 
This case corresponds, in particular, to optical pumping 
by CO2 or quantum-cascade lasers; 

(d) When fifi > huJo (or HV. ^ fiwo), and the character- 
istic times of inter-carrier and optical phonon scattering 
are of the same order of magnitude {taucc ~ tq), the 
pattern of the energy relaxation of photogenerated car- 
riers can be fairly complex. Nevertheless, even in this 
case, on can introduce the energy Sq — hTlo, which, how- 
ever, is determined not only by the ration of hTl and 
hiOQ, but also by the ratio of the characteristic scatter- 
ing times. Naturally, hft — 2Khajo < Eq < hil. In 
this case, the quantity Eq (or fig) is a phenomenologi- 
cal parameter, which determines the fractions of the ab- 
sorbed photon energy going to directly to the electron- 
hole plasma and to the optical phonon system. For a 
rough estimate, the following simple formula might be 
used: H^o = — '2'KhuJo/[l + Ktq/tcc], which yields the 
above values of h^o in the limiting cases (a) and (b). 



Introducing parameter Eq (or J7o), one can consider all 
the abovementioned tracks in the framework of the same 
treatment. 

In the situation under consideration, the optical 
phonon system can also be far from equilibrium, so 
that the distribution function of optical phonons TVo can 
markedly deviated from its equilibrium value A/'q'^ = 
[exp(?iLLio/To) — 1]~^. The latter can lead to an effective 
reabsorption of optical phonons by electrons and holes 
with a significant increase in their energies. Although at 
low temperatures, the energy relaxation and recombina- 
tion of electrons and holes are associated with the inter- 
action with acoustic phonons and radiative processes, re- 
spectively [l^l , in the case of elevated temperatures (say, 
room temperatures) and pumping intensities, the energy 
relaxation and recombination is assumed to be due to the 
interactions with optical phonons (both inraband and in- 
terband) . 

The quasi- Fermi energy of the electron- hole plasma Ep, 
its effective temperature T, and the number of optical 
phonons TVo obey the equations 

hwoiK'"'' + = Gnom, (2) 

fj^^Rdecay ^f^^Q^^ (3) 

governing the balance of the electron-hole pairs, the 
electron-hole plasma energy balance, respectively, as well 
as the optical phonon balance. Here, i?™*'^'', i?™*'"", and 
j^decay -j-j^g ratcs of the pertinent processes. Equa- 
tion (3) explicitly takes into account that all the energy 
received by the system from optical pumping goes even- 
tually to the thermostat (to the contacts via acoustic 
phonons) . The optical generation rate is given by 

Go = 7ratanh( — j I. (4) 

Here, / is the photon flux of pumping radiation and 
a ~ 1/137 is the fine structure constant (so that the 
absorption coefficient is equal to na ~ 0.023). The last 
factor in the right-hand side of Eq. (4) accounts for the 
limitation of the interband absorption associated with 
the Pauli blocking. In this case, Gn — Gjiq — iral. To 
provide an effective cascade pumping and avoid the fast 
direct transitions with the emission of an additional opti- 
cal phonon, the following conditions should be satisfied: 

Ef < < hujo — Ep (5) 

III. RATES OF THE PROCESSES UNDER 
CONSIDERATION 

For the terms i?™*'^'' and R^^^'^ , which describe 
the electron-hole recombination and generation processes 
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and the intraband energy relaxation assisted by optical 
phonons, one can use the following simplified formulas: 



-Dtnter 



nv, 



(TVo + 1) exp 



2ei? — hujQ 
T 



-^^o 



So 



(A/'o + l)exp - 



T 



-Mo 



(6) 



(7) 



Here, Eq = '^(To/fivw)'^ /G is the equilibrium electron 
and hole density and Tq"*'^'' cx Tq"*''" oc tq are the times of 
the interband and intraband phonon-assisted processes. 
The latter quantities have been introduced in the present 
form for convenience 



A difference in r™*'^'' and t, 



intra associated with the 



features of the density of states and with a difference 
between the optical phonon energy and the quasi-Fermi 
energy (see Appendix). In equilibrium, i.e., at A/o = A/'q', 
= 0, and T = To, Eqs. (6) and (7) yield = 
i?™*™ = 0. 

The rate of optical phonons decay due to the anhar- 
monic contributions to the interatomic potential, result- 
ing in the phonon-phonon scattering and in the decay of 
optical phonons into acoustic phonons and is assumed to 
be in the following form: 



is the characteristic photon flux and rj^'^'^"'^ 



decay i inter 



'0 



It™ 



As follows from Eq. (9), at very large 



pumping intensities, TVo can markedly exceed unity. We 
shall not consider the range of such intensities. There- 
fore, in the following we limit our consideration by 
not too strong pumping assuming that 77q"^'^°^///o <^ 

exp(&^o/ro), i.e., / « /o exp(;K^o/To)/r/o'™^ = . 

Using Eqs. (1), (2), (5), and (7) with Eq. (9), we arrive 
at the following equations for the effective temperature T 
and the quasi-Fermi energy ep (for not too strong pump- 
ing intensities): 



f2£p — hajQ 



1 + ivf'' + 1) 
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: exp 
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decay 



(10) 



decay 



So(A/'o 



(8) 



where Tq ^ is the pertinent characteristic time. Con- 
sidering high heat conductivity of GLs ^3], the lattice 
temperature, i.e. the temperature of acoustic phonons, 
is assumed to be equal to the temperature of the contact 
To. 
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, (11) 



IV. CALCULATION OF THE EFFECTIVE 
TEMPERATURE AND QUASI-FERMI ENERGY 



Equations (1), (3), and (8) yield 



(9) 



exp 



l + (%^^''^ + l)^ 
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decay , ( ^0 
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(12) 



Here 



T K req ^0 ^0 

-^0 = 77 i;it57 - exp 



To J n nar^" 



Here 770 = r^"t™/r™*«'-. 

Equations (11) and (12) yield 
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l--^lnO 
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decay 

Vo 



(^0 , 

Vo\ 1 

Vwo 



(13) 



As it should be, Eqs. (13) and (14) yield T = To and 
ep, = if / = 0. 

The dependence of parameter 779 on sp and T compli- 
cates the solution of Eqs. (13) and (14) in wide ranges 
of pumping intensities, particularly, their analytical so- 
lution. 



V. SPECIAL AND LIMITING CASES 

A. Special case fio/'^o = !• 

In a special case ilo/^o = 1 (i-C-, ^ = {^K + l)a;o), 
Eqs. (12) - (14) yield 



T : 



To 



l_i<L ln(l+^f-^///o) 



n 



To + ^ ln(l + 77f™^///o), (15) 
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In 1 
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-^Ind 



decay t I t \ 



I/Io 



-1 I decay t / t 

1 + "I/Io 



exp 



2eF 



= 1 



I/Io 



(16) 



(17) 



Thus, at fio/wo = 1, the effective temperature T is an 
increasing (logarithmic) function of the pumping inten- 
sity /, although the increase in T is rather slow due to a 
small factor To/fiwo — 0.125 in the last pre-logarithmic 
term in the right-hand side of Eq. (15). As follows 
from Eq. (16), the quasi-Fermi energy ep logarithmi- 
cally increases with / and tends to the saturates with 
= {To/2) ln(l + 1/r]^'""'^). The increase in e^. in the 
range of moderate I/Iq (I/Io ^ ^/Vdecay) is faster than 



that in T. One can also see that an increase in the decay 
time leads to slowing down the sp — I dependence. 

One can see that the variation of the effective tem- 
perature (the electron-hole plasma heating) is associate 
with the finiteness of the optical phonon decay time, i.e., 
with the optical phonon heating and the absorption of 
hot phonons by electrons and holes, while the contribu- 
tions of the intraband and interband emission and ab- 
sorption of optical phonons to the electron-hole plasma 
energy balance at i^o/wo are compensated. It is instruc- 
tive that in the case in question, the effective tempera- 
ture of the optical phonon system O, defined in such a 
way that TVo — [exp(fia;o/6) — exceeds the lattice 
temperature To (optical phonon heating) and is equal to 
the effective temperature of the electron-hole plasma T. 
Indeed, using Eq. (9), we find (compare with Eq. (15)) 



e 



To 



1 - 



tvjJo 



ln(l 



Vt'^'l/h) 



(18) 



The fact that 8 = T > Tq implies that both electron- 
hole plasma and optical phonon system are heated and 
they are in equilibrium with each other. 

If Ml — 1 eV, so that K — 2, assuming that To = 
300 K, and rj^'^'^"'^ = 1 — 3, and setting I/Iq = 1 / tuq'^'^"'^ , 
we obtain T 9 ~ 326K, When I/h > 1/770"™^, we 



(7 - 17) meV, and exp(£|?72r) < 



obtain ep — £p 
1.07-1.19. 

The optical phonon decay time and, hence, parameter 
^decay j^jg]^^ ]-,g gj^^U in GLs On properly chosen sub- 
strate. In such a case, the quasi-Fermi energy at suf- 
ficiently strong pumping can be not so small. Indeed, 



setting I/Io = ^/v, 



decay 




and 



^70 



ecay 
< 



= 0.1 — 0.5, from 
17 — 45 meV, and 



Eqs. (15) - (17) we obtain ef^ 
exp(e|?72r) < 1.19 - 1.57. 

The optical phonon decay time, particularly in sus- 
pended GLs, can be fairly long (in the range 1 — 10 ps 
[23I . [31I - I33I I ) , and it might markedly exceed the character- 
istic time of intraband interaction of electrons and holes 
with optical phonons. Hence, the inequality r\^'^'^^ ^ 
ryo > 1 can, be valid, particularly, in the case of sus- 
pended GLs. In this fairly realistic case, Eqs. (12) - (14) 
lead to Eqs. (14) - (17) and their consequences at all val- 
ues of rio/wo, because the terms with 770(^10/^^0 — 1) are 
relatively small and can be omitted. 

Thus, at long optical phonon decay times, the effect 
of accumulation of optical phonons (heating of the opti- 
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cal phonon system) can prevent realization of population 
inversion necessary for lasing (large quasi-Fermi energy 
ep) at room temperatures when the processes involving 
optical phonons can dominate. 



B. Weak pumping 

At sufficiently weak pumping (see below), Eqs. (12) 
- (14) can be solved analytically at arbitrary values of 
?7g'°^"^, rjo, and fio/wo- In this case, the effective temper- 
ature and quasi-Fermi energy are close to To and zero, 
respectively. As shown in the Appendix, a,t ep <T, one 
obtains r]Q = ~ (fiwo/7rro)2/(l + 2.19To/huJo). At 
room temperature, jjg^ ~ 5. Thus, at a weak pumping. 
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Simultaneously, one obtains 



exp(^).l + 



no 

Wo 



- 1 



To 



6 Wo / lo 



(19) 



(20) 



(21) 



However, at sufficiently small ratio f2o/wo, parame- 
ter ■q'^'^ can be negative, and the effective temperature 
can decrease (cooling of the electron- hole plasma). If 
f2o/wo < 6/5, the quasi-Fermi energy is a rising function 
of the pumping intensity. The latter inequality implies 
that there should be $1 < uJo{2K + 1.2). The latter con- 
ditions of the electron-hole plasma cooling accompanied 
with the quasi-Fermi energy increase do not contradict 
inequalities Eq. (5). The effect of cooling under consid- 
eration is attributed to the effective photogeneration of 
low-energy electrons and holes (with the energy f£lo/2) 



decay 



C. Strong pumping, fio/wo > 1 — Vo' ho 

In this case, ■q'^'^''"'^ + r]^{^lo/LOo — 1) > 0, even at el- 
evated pumping intensities / > /q, a moderate increase 
in T and exhibit a moderate increase with increas- 
ing I/I01 which slows down in the range of large I/Io 
(the T — I dependence becomes logarithmic, whereas the 
quasi-Fermi energy tends to a constant value, positive 
or negative). Indeed, in the case under consideration, 
Eq. (14) does not have solutions with ep '^T. Hence in 
the logarithm in Eq. (14) one can put % = fff — 5. As 
a result, we obtain 



T~Tf 







rp2 



In 



To+f^lln 
\njJo 



^0 1 
Ew — — In 

2 



decay 



Wq 



r^-1 



\Wo 
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Equations (19) - (21) are valid if 



^decay_^_^eql^_-^ 
Wo 



\V' 



where rj""^ = vi},'""''" + Vjl'' {^lo / ujq - 1). 

As clearly seen from Eqs. (19) and (20) (as well as 
from Eqs. (13) and (14)). the effective temperature in- 
creases (the electron and hole heating) with increasing / 

if r?^9 > 0, i.e., if f2o/wo > 1 - %"'"7'?o''- The rise of 
the effective temperature is more pronounced at larger 
^decay ^ i.e., at longcr decay time Tq^'^"'^ . Under the con- 
dition ilo/wo > 1 + l/??o^ ~ 6/5, the quasi-Fermi energy 
can become negative ep < 0, so that cxp{ep/T) < 1. In 
the latter case, the optical pumping does not lead to the 
degeneration of the electron-hole plasma. It is interesting 
that this conclusion is independent of the relative value 



of the decay time Tq 



decay 
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^decay ^ ^ 



decay 
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When 77^^"^ » 1, Eq. (23) yields 



.sat ^ ^ In 



1 + 



^decay y 5 



^0 
Wq 



2% 



no 

Wq 



(23) 



(24) 



As follows from Eqs. (20), (23), and (24), ep changes its 
sign at fio/wo =6/5 both at weak and relatively strong 
pumping. In the case under consideration, the popula- 
tion inversion is weak {e p can be positive , but small) or 
is not realized {sp < 0). 
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Formally, the logarithmic factors in Eqs. (22) and (23) 
diverge when ^Iq/cjq tends to 1 — tJq'^ / rj^^''"'^ ■ However, 
these equations become invalid in the immediate vicinity 
of the point fJo/'^o = 1 ~ '7o^/'?o^'^°^: because above we 
have assumed ep/T < 1 and put r]o = t^q'. The case of 
smaller flo/^o is considered in the next subsection. 



at sufficiently large I /Iq. In such an instance, one can ex- 
pect that ep ^ T and r/o — {hujo/ep)/6. Hence, Eq. (14) 
can be presented as 



D. Strong pumping, small fio/wo (cooling regime). 

The case ^Io/loq is small, is much more interesting be- 
cause rather large values of ef and ep/T can be achieved 
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(26) 



r 



Assuming that ep ^ T, one can approximately solve 
Eqs. (25) and (26) by iterations, considering the small- 
ness of exp(— £ir/T), and arrive at 



T~To 



1-2 



Wo / /o 



6 l + ??o 



(27) 



\ 



Wo / lo 



61 l + rj^^^'^y- 



(28) 



If I/Io > 1/%'"'"^. from Eqs. (27) and (28) 
for T and ep we obtain their satu ration values: 

T-* ~ To[l - 2^(l-Q o/c.o)/6r;o^^n and ef' ^ 

TujJq^J (1 — Oo/wo)/6?7g"^"^ . The corrections to Eqs. (27) 
and (28) are exponentially small when is not too 

close to unity. Due to the above assumptions, Eqs. (27) 
and (28) are valid until 2eF < hoJo, that yields 



' 2To 



2(1 - no/ujo)I/Io 



3(1 + %^^^"^)///o 



< 1. 



Considering that (2To/3;ia;o)^ - 4/576 < 1, these equa- 
tions are valid in rather wide range of pumping intensi- 
ties, (not too small), except the case when (1 — fJo/wo) 
is very close to zero (this is the special case considered 
above). Some complications can arise in a non- realistic 
case of very small t^q'^^"*'. Using Eqs. (27) and (28), one 
can find the condition when ep ^ T that was assumed in 
deriving Eqs. (25) and (26) and their consequences (see 
also Eq. (5)): 



^<1 

Wo 



a 

6% 



decay 



\hujQj 



At room temperature, the latter inequality yields 
fio/wo < 1 — O.ljyg'^''''^. This implies that large ep and 
T < To can be achieved if tiq^'"''^ < 10 and the ratio 
fio /wq is chosen properly. 

As it was mentioned above, the effect of cooling under 
consideration is attributed to the effective photogenera- 
tion of low-energy electrons and holes. 

Setting To = 300 K, Ml = 900 meV, we obtain 
HQq — 0.1 eV and f2o/wo — 0.5. For this case, assum- 
ing r]^'""'^ = 1 - 3, we obtain T"»* 127 - 200 K and 
ef* ~ (33 - 58) meV. Hence, exp(e|,»72r) ~ 2.7-15.5. 
These latter two values are markedly larger than those 
obtained in Subsection A for fio/wg = 1. If ?7q'^'^"^ < 1, 
the effect of cooling is more pronounced and the quan- 
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FIG. 1: Normalized carrier eflective temperature T/Tq versus 
normalized pumping intensity I/Io for different r/g'^'^''^ and 
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FIG. 2: Quasi- Fermi energy sf versus normalized pumping 



intensity ///o for different r\^'^°'^ and Ho /wo 



titles e|?* and exp(e|?*/2r) can be fairly large. At the 
parameters chosen for the above estimates, conditions (5) 
are satisfied. 

It is instructive that In the cooling regime in question, 
the optical phonon system is heated. The effective optical 
phonon temperature in this regime is given by the same 
equation as in other regimes under consideration, i.e. by 
Eq. (18), which yields 8 > Tq. 



VI. CALCULATION OF THE GL DYNAMIC 
CONDUCTIVITY 

The attenuation and amplification of electromagnetic 
waves and surface plasmons propagating along the GL 
plane is determined by the the real part of the GL dy- 



namic ac conductivity at the signal frequency w (photon 
or plasmon frequency) Re tXi^. The latter includes the 
contributions of the interband and intraband transitions: 



Re a^i — Recrf 



Re erf 



(29) 



For the multiple-GL structures (see, for instance, Refs. 
(U), the contributions of all GL to Re should be sum- 
marized [U, [l2| . 

In the situation under consideration, generalizing the 
pertinent formula from Ref. @ (see, for instance, Refs. 0, 
[ll|) for the photogenerated nonequilibrium electron hole 
plasma with the equal non-zero quasi-Fermi energies of 
electrons and holes, one can obtain: 



Rear 



{^— \ tanhf 

V4^y V 



4T 



(30) 



Here, e— |e| is the electron charge. At w ^ r^^, where 
T is the momentum relaxation time. The intraband con- 
tribution, which actually corresponds to the Drude ab- 
sorption by electrons and holes, can be obtained from the 
Boltzmann equation 0, |Tl|, [s^ : 



RecT,* 



2e2 



dee d 
T de 



cxp[(£ - £f)/T] + 1 

(31) 

As To and T are small in comparison with hujo , the elec- 
tron and hole momentum relaxation is associated with 
different scattering mechanisms, in particular, with the 
scattering on short-range and long-range disorder and 
acoustic phonons. In this case, = voie/To), {vq is the 
collision frequency of electron and holes with the thermal 
equilibrium energy), so that interpolating the integral in 
Eq. (31) [see Eq. (AlO)] as a function of e and T, from 
Eq. (31) we obtain 



Recr, 



intra 



2e^iyo (el + n^Tye) 



To 



When > T from Eqs. (31) and (32) we obtain 

T3„ intra ^ 2e^VQ e\ 
^^^^ - fc2 

Hence, using Eqs. (30) - (33), we arrive at 



(32) 



(33) 



Re CTto = — 



/ hoj — 2eF 

tanh 

V 4T 



(4 + 7r2TV6)- 



(34)" 

One can see that at r ^ oc e, the Drude term given by 
Eqs. (34) and (35) is proportional to the density of the 
electron-hole plasma. 

In the special case considered above when rio/ujo = 1, 
using Eqs. (16), (17), and (34), we obtain 



8 



Re (T[, 



exp(/i(i;/2ro) 



(1 + ryf ""'«///o)-/2"o 



- 1 



1 . decay j i j 

1 + ?7o ^///o 



1/2 



exp(/ki;/2ro) 



(1 + r;^f"™^///o)"/2"" 



1 . decay J I J 
1 + % ''///o 



1/2 



w / 1 27r2 



1 I decay t I t 

1 + n/h 



(35) 



Here ctq = e^/4?i and w_d = ^47rt'oTo/3?i. In the THz and middle IR ranges, a; < 2wo (w/27r <C 10^ THz), Eq. (35) 
can be reduced to the following: 



i/h 



1 I decay j i j 

1 + nji^ 



1/2 



1 I decay j i j 

1 + Ujl^ 



1/2 



decay t i t 



r 



1 - ln(l 



decay 7- / 7- 



(36) 



Assuming j/o = (0.25-1.0) x lO^^ g-i Tq = 300 K, we 
obtain fo = a;£)/27r ~ (1 — 2) THz. Considering that the 
momentum relaxation time in multiple-GLs at the tem- 
peratures To < 50 K can reach values about 20 ps [33 |. 
the above values of vq at Tq = 300 K appears to be rea- 
sonable (assuming a linear increase of with the lattice 
temperature Tq [33 )• 

When f^o/wo is smah, using Eqs. (27), (28), and (34), 
we arrive at 



Re (Ju_ 



exp 


hjj 1 


— exp 


hwQ Z 


[2ro (1-2Z)J 


[ To {1-2Z)\ 


exp 


hijj 1 


+ exp 


hijjQ z 


_2Tq (1-2Z). 


[ To (1-2Z)J 





2 








) 


_7r2 


kTo J 



(37) 



Here Z ^ J [{1 - ^o / luo) I / h] / {1 + vfl/Io) 



VII. DISCUSSION OF THE RESULTS. 

As follows from Sec. V, to achieve the population inver- 
sion sufficiently strong for lasing, the parameter Hq/luo 
should be either equal to unity (the special case) or 
smaller than unity (cooling regime) . Below we shall focus 
on this situation. 



Figure 1 shows the dependences of the normalized ef- 
fective temperature T/To on the normalized pumping in- 
tensity I/Iq calculated for different 1]^'"^°'^ using eq. (15) 
for ^lo/<^Q = 1 and Eq. (27) for fto/coo — 0-5. Figure 2 
shows the dependences of the quasi-Fermi energy sp on 
the normalized pumping intensity I / lo calculated for dif- 
ferent tjq'^'^"'^ and ^lo/^o using Eqs. (16) and (28). 

One can see from Figs. 1 and 2 that at ilo/uto < 1, 
a pronounced cooling of the electron-hole plasma takes 
place, and this cooling regime provides larger sp for a 
given pumping intensity. As also follows from Figs 1 
and 2, the finiteness of the optical decay time markedly 
affects the T — I and ep — I dependences. Indeed, at 
r^o/wo — 1 , the T — I dependence becomes steeper when 
the parameter j]q'^'^"'^ increases (reinforcement of heat- 
ing effect), whereas at Qq/ujo < 1, this dependence be- 
comes less pronounced (weakening of cooling). This is at- 
tributed to the optical phonon heating, which promotes 
the heating of the electron- hole plasma at Qq /ujo = 1 and 
counteract its cooling at flo/too < 1- 

Figure 3 demonstrates the dependences of the real part 
of the GL dynamic conductivity Re (Ti^ on the signal fre- 
quency ll)/2tt at different normalized pumping intensities 
I/Iq calculated for different rj^'^'^"'^ and vo using Eq. (35) 
(assuming that f2o/wo = 1). As seen from Fig. 3, at 
rio/ti^o — li the dynamic conductivity Re cr^j can be neg- 
ative only if parameter rj^^'^"'^ 1 (the optical phonon 
decay time is small), i.e., if the heating of electrons, holes, 
and optical phonons is very weak. However, at VLo/ojo < 1 
when electron-hole plasma is cooled. Re (Jt^ < in range 
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FIG. 3: Real part of the normalized dynamic conductivity as a 
function of signal frequency at different normalized pumping 
intensities I/Io and different t^q'^'^"^ ( Qo/^o ~ 1 and uo = 
0.25 X fO^^ s"^). 



of signal frequencies from lu/2tt ~ 2 THz to 17 THz and 
more [see Figs. 4 and 5, obtained using Eq. (37)]. This 
occurs at not too low pumping intensities (I/Iq — 5) 
but even at rather large values of i'q and tjq^'^'^^ when the 
cooling of the electron- hole plasma is accompanied by the 
heating of the optical phonon system. 

Figures 6 and 7 calculated using Eq. (37) show the 
real part of the GL dynamic conductivity Re a^j ver- 
sus the normalized pumping intensity I/Iq calculated 
for given values of the signal frequency uj/27t — 2 THz 
(at J/o = 2.5 X 10" s-i) and w/27r = 3.5 THz. (at 
1^0 = 1 X 10^^ s~^). In both cases, Re is negative 
when I/Iq exceeds certain threshold value. The latter is 
somewhat smaller in GL with smaller vq . It is interesting 
that in the cases v^'^'^"^^ = 1, an increase in I /Iq leads to 
smaller |Re (Jt^\ in the range of large ///q, in contrast 
to the cases z/q^'^"*' = 2 and v^'^'^"'^ = 3, This can be ex- 
plained by a faster increase in ep with increasing ///q at 
smaller Vq'^'^'^^ leading to a faster increase in the Drude 
absorption. 

To estimate Iq, we use the data for Sq/tq"*'^'' at 
To — 300 K from Ref. [lB| (see Appendix) and set 
hQo = 100 meV, so that (fio/wo — 0.5). and setting 
hrio = 100 meV. The chosen value of hflo can correspond 
to the optical pumping photon energy hQ = 900 mcV 
provided the cascade emission oi K — 2 optical phonons 
or to h^l = H^Iq = 100 meV (pumping by CO2 laser). 
In such cases, one obtains Iq ~ 9,66 x 10^^ 1/cm^s and 
Iq ~ 8.69 X 10^^ 1/cm^s, respectively. In both cases, 
^hfllo- 1.39 X 10* W/cm2. 

As demonstrated, lowering of the ration Qq/ujo is bene- 
ficial for achieving the population inversion and negative 
dynamic conductivity, because such lowering results in 
weakening of the electron-hole plasma heating and even 
in the transition to its cooling. Moreover, as follows from 
Eq. (23) and (24), at elevated flo/ujQ, the value of the 



FIG. 4: Real part of the normalized dynamic conductivity as 
a function of signal frequency at different normalized pumping 
intensities I/Io and different tJq'"^"'^ ( Qq/ljo = 0.5 and uo = 
2.5 X 10" s"^). 
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FIG. 5: The same as in Fig. 4 but for = 1.0 x 10^^ s"^ 



quasi-Fermi energy ep can become negative that corre- 
sponds to nondegenerate electron-hole plasma (/ < 0.5). 
In this regard, in the case of optical pumping with rel- 
atively high photon energy, the suppression of the pair 
carrier interactions might be indispensable for the real- 
ization of strong population inversion and lasing. This, 
possibly, can be done using the GL structures with high-k 
substrates. 

Instead of the optical pumping under consideration, 
the electron-hole cooling and the realization of strong 
population inversion can be achieved by the injection 
pumping (in GL structures with p-n or p-i-n junctions) 
considered previously Q. In such a case, one can put 
rio — eV/h <C Wo, if the applied bias voltage V is not 
large {V < 200mV). 

Above we assumed for simplicity that there is only one 
type of optical phonons with Hujq ~ 200 meV. Actually, 
the processes involving the optical phonons with uji = 
loq) and uj2 = ujq — Awq are important ill (Awo/wQ ~ 
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FIG. 6: Real part of the normalized dynamic conductivity as 
a function of normalized pumping intensities //Jo for lojl-K — 
2 THz and different r]^^'^'^^ ( Q.o/ujo = 0.5 and i^o — 2.5 x 
10^1 s-i). 



Eq. (34)]. Therefore, although, as demonstrated above, 
the THz and IR lasing using the population inversion 
in graphene can be realized even at room temperatures, 
at lower temperatures, it can be achieved easier. How- 
ever, at lower temperatures, the optical phonon assisted 
processes in GLs can give way to other recombination, 
generation, and energy relaxation mechanisms, for in- 
stance, the radiative recombination and generation and 
the energy relaxation on acoustic phonons (l3l [l3 . |36| . To 
estimate the temperature range, in which the above sim- 
plified "optical phonon" model is sufficient, we assume 
that the characteristic time of radiative recombination 
and that associated with the optical phonon emission are 
about (10 - 20) ns and 0.5 ps, respectively. In this case, 
one can obtain To > 220 — 230 K, so that our model is 
valid at room temperatures and slightly lower. 



VIII. CONCLUSIONS 
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FIG. 7: The same as in Fig. 6 but for for cj/27r = 3.5 THz 
and I/O = 1.0 x 10^^ s'^ 

0.18). This adds complexity to the pattern of the electron 
and hole relaxation affecting the value of parameter f^o 
and 

Although one can expect that due to a smallness of 
Awo/'^Oj this effect is not essential and can be accounted 
for by a proper correction of fio- • More complicat- 
ing factor can be the interaction of electrons and holes 
with optical phonons in the substrate. These phonons, 
having relatively low energy, can, in principle, markedly 
contribute to the recombination and energy relaxation. 
Perhaps in multiple GL-structures the effect substrate 
phonons is not pronounced. The pertinent generalization 
of our treatment is going to be considered elsewhere. 

Lowering of the lattice temperature Tq and, conse- 
quently, the effective temperature T results a weaken- 
ing of the efficiency of optical phonon assisted processes, 
particularly, in a decrease in the recombination rate [see 
Eqs. (6) and (7)]. Apart from this, the Drude absorp- 
tion also becomes smaller due to a decrease in vq [see 



We have considered the characteristics of electron-hole 
plasma in optically pumped GL at elevated (room) tem- 
peratures (its quasi-Fermi energy, effective temperature, 
and dynamic conductivity) and shown that: 

(1) The interband and intraband processes of emission 
and absorption of optical phonons play a crucial role 
in the characteristics of optically pumped electron-hole 
plasma in graphene at elevated temperatures; 

(2) The electron-hole plasma in optically pumped 
graphene can be both heated and cooled depending on 
parameter fig/i^oj while the optical phonon system is al- 
ways heated; 

(3) The accumulation of nonequilibrium optical phonons 
(their heating) due to finiteness of their decay time leads 
to a slower dependence of the electron and hole quasi- 
Fermi energy on the optical pumping intensity; 

(4) The effects in question can markedly influence the 
achievement of the negative dynamic conductivity in op- 
tically pumped GLs associated with the population inver- 
sion and, hence, the realization of THz or IR lasing in the 
GL structures at room temperatures. The latter requires 
a careful choice of parameter fJo /i^o and minimization of 

, decay 

parameter r^Q 



IX. ACKNOWLEDGMENT 

The work was supported by the Japan Science and 
Technology Agency, CREST, Japan. 



Appendix 

For the term i?™*'^'' in Eq. (6), which describes the 
electron-hole recombination and generation processes as- 
sisted by optical phonons, one can use the following sim- 
plified formula: 
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Jointer 



T3 , 



hujQ 



dee{hjjQ — e) 



1 + exp 



e — ep 
T 



1 + expl 



A/-o + l-exp(^^^^^L^)A/-o 



6T3 



(M, + l)exp(^£^::^)-A/-o 



du 



oxp(-eF/T) 



(l + u)[l + ue 



(eF-huJo)/T^ 



6T§ 



{Mo + 1) exp 



2£p — hujQ 
T 



( hujQ — 2£f 

exp I 



T 



exp 



T 



In^exp ^- 



T 



exp\ I + l-,2 



(Al) 



Equation (Al) can be markedly simplified in the most 
interesting situations when htoo > 2ei?,T: 



with 



(A/-o + l)exp(^£^^)-AAo 



So 



(A2) 



(A3) 



Here we have neglected the terms in the pre- 
exponential factor with the ration Sp/hujQ and in- 
troduced the rate of the thermal generation of the 



electron-hole pairs due to absorption of optical phonons 
in equilibrium Gg"^ = Go{hujQ)^ /6Tq exp{—hujo)/T) = 
(So/r™*"^'') exp(— /iwo/T)- The latter was estimated in 
Ref. 161 



10^^ cm-^s"^ at To = 300K. 



Equation (A2) differs from those in Refs. [ll|, [l^ by 
the inclusion of the terms with A/q, which correspond to 
the processes of stimulated emission and absorption of 
optical phonons At temperatures, T ^ hojQ ~ 2300 K, 
when the number of equilibrium optical phonons A/q = 
A/'o'' <^ 1 (stimulated emission and absorption of optical 
phonons^is negligible), eq. (A2) reduces to that obtained 
in ref. 

The term i?™*''" in Eq. (7) can be calculated as 



Rn 



2Go 



dee{e 


+ huja) 




+ 1 exp ^ 




So 


1 + exp 1^ 


' e — ep 


)] 




1 -1- exp 1 




—intra 


T 






T 







(A/'o + 1) exp 



(A4) 



where 

So 



1 + exp(-it ~sp/T) 



If IefI < T, Eq. (A4) yields 



So 



TT^Go f tY fhujo + 2.19T\ 



6 VTo 



To 



(A6) 



Comparing Eqs. (A3) and (A6), for this case we obtain 

.o=^=4r^Vfi+2.i9-^y\ (A7) 



If ftwo > > T, 

So ^ fujoel,{l + 2eF/3fmjo) fujoej, 

^0 ^0 ^0 

so that invoking Eq. (A3), we obtain 

2 



1 ( hujQ 



6 V ef 



(A9) 



In a wide range of variations of Ef and T, one can use 
the following interpolation: 



Vo = 



(64 + 7r2T2)' 



(AlO) 
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